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if  the  events  at  the  two  sites  are  perfectly 
correlated  and  mutually  dependent. 

If  we  denote  the  absence  of  the_ event  at 
site  A as  P(W^),  and  at  site  B as  P(WB),  then 
P (WA)  =1-P  (H^)  and  P(WB)ml-p(^)  , and  we  can  write 
Eq.  1 as: 

P(WaWb)=  [l-P  {vTA7j  [l-P  (W^TJ 

= l-Pj^-PfW^+PtW^PtW^)  (3) 

Eq.  3 applies  only  in  the  case  when  the 
events  are  statistically  independent.  For  the 
case  when  the  events  are  not  statistically 
independent,  the  joint  probability  can  be  esti- 
mated with  the  expression 

P (WAWB)  =1-P  (WA>  -P  (Wb)  +P  (»Ta)  P (wb>  Kl  (4) 

P(WA)>  P(WB) 


INTRODUCTION 


Meteorological  satellites  have  made  it 
possible  to  obtain  information  on  areal  coverage 
and  joint  occurrences  of  meteorological  events 
not  heretofore  possible.  However,  a long  record 
of  observations  is  required  to  derive  a climatol- 
ogy of  joint  occurrences.  Routine  hourly  surface 
weather  observations  may  be  sufficiently  con- 
current to  estimate  joint  occurrences  among 
sites  where  such  observations  are  taken.  Satel- 
lite observations  are  especially  useful  in  areas 
where  surface  observations  are  unavailable  or 
unsuitable . 


Unconditional  probabilities  of  weather 
events,  such  as  precipitation,  freezing  temper- 
atures, overcast  skies,  etc.  can  be  estimated 
with  considerable  accuracy  from  available  records 
However,  the  probabilities  of  such  events  occur- 
ring jointly  at  two  or  more  locations  is  often 
difficult  to  estimate  if  the  events  are  not 
statistically  independent.  Joint  probabilities 
of  a weather  event  can  be  estimated  either  from 
relative  frequencies  obtained  directly  from  the 
data  or  from  a model.  Large  samples  of  data  must 
be  processed  to  obtain  relative  frequencies  that 
are  good  estimates  of  true  probabilities.  Models 
eliminate  the  need  for  data  processing  to  esti- 
mate joint  probabilities  and,  in  some  cases,  they 
can  actually  provide  better  estimates  of  true 
probabilities  than  those  obtainable  from  data 
samples,  except  in  those  cases  where  the  data 
samples  are  very  large. 


When  there  is  perfect  correlation,  Ki“0  and 
Eq.  (4)  reduces  to  Eq.  (2) . When  the  events  are 
statistically  independent,  Kx”l  and  Eq.  (4) 
reduces  to  Eq.  (3).  Nearly  perfect  correlation 
might  be  expected  when  the  sites  are  very  close 
to  one  another.  Statistical  independence  is 
approached  as  the  distance  between  the  sites 
becomes  large.  Obviously  is  a function  of  the 
distance  between  sites  A and  B. 


2.2  Three-site  Joint  Probabilities 


When  weather  events  are  statistically 
independent,  three  site  joint  probabilities  can 
be  expressed  as  follows: 


P(WAW'  [l-P <"AJ  [l-P  0*B>J  [i-PWc>J  <! 

When  the  events  are  not  independent  they  can  be 
estimated  with  the  following  expression: 


2.1  TWo-slte  Joint  Probabilities 


Let  us  denote  the  probability  that  a 
weather  event,  W,  will  occur  at  site  A by  P (WA) , 
at  site  B by  P(WB),  and  jointly  at  the  two  sites 
by  P(WAWB) . We  know  that 


P <mAnBnC> "l-P <WA) "P (*B> "P  <«C>  +P  <*A> p <*B> 
+P  <WA)  P (WC)  Kl+P  (Wp)  P (tTc) Kl 
-p(wA)P(w^)Kip<wc)K2 

P(ffA)*  P(5b)  £ P(ifc) 


P <"aWB  ) *P  <WA)  p <»B ) U) 

if  the  events  at  the  two  sites  are  statistically 
independent . 


We  also  know  that 


The  ICj  values  are  obtained  from  the  same 
function  as  those  used  in  Eq.  (4) . The  distant 
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between  sites  A and  B,  AB,  is  the  only  distance 
required  for  finding  Kj  in  Eq.  (4) . Distances 
HE,  Sc , and  Sc  are  required  for  determining  the 
necessary  X^  values  in  Eq.  (6)  and  Kj  is  a func- 
tion of  the  distance  SC  or  SC,  whichever  is 
shorter . 


2.3  N-site  Joint  Probabilities 


In  the  general  case,  say  n sites,  joint 
probabilities  of  independent  events  can  be 
expressed  as  follows: 

P (WAWB . . .W  > = li-P  (wj  fl-PlwJ...  fl-P  rw_n 


When  the  events  are  not  independent,  they 
can  be  estimated  with  the  following  expression: 


Figure  1.  Location  of  the  nine  stations  whose 
winter  and  summer  hourly  observations  of 
precipitation  were  studied. 


In  winter,  the  observers  recorded  0.0641, 
0.636  and  12.638  percent  of  the  observations  as 
heavy,  heavy  or  moderate , and  all  intensities  of 
precipitation,  respectively.  In  svsgmer,  the  per- 
centages were  0.202,  0.4S9  and  6.257.  Clearly 
most  of  the  precipitation  events  were  of  the  light 
or  very  light  variety. 


is  found  for  the  distance  between  the 
two  sites  with  the  highest  probability  of  the 
weather  event  not  occurring,  in  each  term  of 
Eq.  (8)  in  which  appears;  K,  is  found  for  the 
distance  between  the  site  witlithe  third  highest 
probability  and  the  closer  of  the  two  sites  with 
higher  probabilities  of  the  weather  event  not 
occurring,  and,  so  on.  The  last  K in  Eq.  (8), 
Kg,,  is  found  for  the  distance  between  site  n, 
the  site  with  the  lowest  probability,  and  its 
nearest  neighbor. 


This  study  addresses  neither  the  subject  of 
precipitation  intensity  nor  the  amount  of  precipi- 
tation. The  joint  occurrences  studied  are 
occurrences  of  active  precipitation  at  two,  or 
more , locations . 


3.1  Unconditional  Probabilities 


Records  of  hourly  precipitation  occurrences, 
observed  in  winter  and  summer,  during  the  13  year 
period  1951  through  1S»3,  at  the  following  nine 
observing  sites,  shown  on  Fig.  1,  were  studied: 


The  relative  frequencies  of  precipitation 
as  reported  on  the  hourly  observations  were  used 
as  estimates  of  the  unconditional  probabilities. 
Because  these  relative  frequencies  were  based  on 
28,060  hourly  observations  taken  at  each  station 
in  winter  and  28,704  hourly  observations  in  siamaer 
they  are  believed  to  be  good  estimates  of  the  true 
probabilities.  The  relative  frequencies  are  shown 
in  Table  1. 


LGA-LaGuardia  Airport,  New  York,  NY 
JFK-Kennedy  International  Airport,  New  York,  NY 
EWR-Newark  Airport,  NJ 

PHL-Philadelphia  International  Airport,  PA 
BAL-Baltimore-Washington  International  Airport,  MD 
DCA-National  Airport,  Washington,  DC 
ADW- Andrews  AFB,  MD 
RIC-Byrd  Field,  Richmond,  VA 
RDU-Raleigh-Durham  Airport,  NC 


Table  1.  Relative  frequencies  of  precipitation, 
RF (W) , at  each  of  the  nine  sites  under 
study. 


Each  hour,  approximately  on  the  hour,  a 
weather  observer  at  each  of  the  above  sites  went 
outdoors  to  make  his  regular  hourly  observation. 
One  of  the  weather  elements  that  he  recorded  was 
precipitation.  The  Federal  Meteorological  Hand- 
book (1975)  describes  how  the  observations  are 
taken. 


For  this  study,  only  active  precipitation 
occurring  at  the  time  of  observation  was  con- 
sidered as  a precipitation  event.  All  types  of 
precipitation  (rain,  snow,  sleet,  hail,  etc.)  and 
all  intensities  (very  light,  light,  moderate  and 
heavy)  were  included. 


SITE 

WINTER 

SUMMER 

LGA 

.1470 

.07643 

JFK 

.1405 

.07639 

EWR 

.1494 

.07624 

PHL 

.1379 

.07238 

BAL 

.1353 

.06881 

DCA 

.1289 

.06730 

ADW 

.1259 

.06487 

RIC 

.1297 

.06951 

RDU 

.1165 

.06970 

3.2  Two-site  Joint  Relative  Frequencies 

As  previously  stated  the  probability  of 
precipitation  occurring  concurrently  at  two 
observing  sites  is  a function  of  the  distance 
between  sites.  The  distances  between  sites  are 
shown  in  Table  2. 


Table  2.  Distances  between  sites,  in  nautical 
miles. 


JFK  EWR 

PHL 

BAL 

DCA 

ADW 

RIC 

RDU 

L<JA  9 

16 

95 

185 

210 

214 

292 

431 

JFK 

21 

94 

184 

208 

213 

289 

427 

EWR 

80 

169 

195 

198 

279 

417 

PHL 

90 

115 

119 

199 

337 

BAL 

29 

30 

121 

256 

DCA 

9 

93 

228 

ADW 

95 

227 

RIC 

138 

The 

relative  frequencies  of  joint  occur- 

rences  of  precipitation  at 

the  36  pairs  of 

observing  sites 

are  shown 

in  Table  3. 

Because 

the  data  samples  are 

large  the  relative  frequen- 

cies  should  be  good  estimates  of  true 

probabilities. 

Table  3. 

Relative  frequencies  of  joint  occur- 

rences  of  precipitation,  RF 

(MW) 

, at 

each  of  the 

thirty-six  pairs  of 

obser- 

vation  sites,  in 

wint< 

ir.  Estimated 

joint  probabilities,  f 

?(WW)  , 

obtained 

from  the  model  are  shown  in 

parenthesis 

• 

J r% 

m 

ML 

ML 

MW 

DC* 

UC 

MJU 

LG*  .125 

.128 

.0942 

.0738 

.0843 

.0847 

.0551 

.0380 

(.128) 

(.124) 

(.0905) 

(.0710) 

(.0859) 

(.0855) 

(.0542) 

(.0385) 

.119 

.0925 

.0718 

.0825 

.0827 

.0538 

.0348 

(.121) 

(.0889) 

(.0881) 

(.0834) 

(.0829) 

(.0523) 

(.0353) 

** 

.0991 

.0771 

.0870 

.0878 

.0579 

.0378 

(.0982) 

(.0751) 

(.0895) 

(.0892) 

(.0588) 

(.0385) 

ML 

.0921 

.0808 

.0808 

.0889 

.0425 

(.0882) 

(.0795) 

(.0788) 

(.0840) 

(.0442) 

.104 

.108 

.0793 

.0518 

(.105) 

(.105) 

(.0771) 

(.0534) 

.108 

.0821 

.0528 

(.112) 

(.0800) 

(.0535) 

.0828 

.0539 

(.0798) 

(.0548) 

.0898 

(.0897) 

3 . 3 Two-site  Estimated  Joint  Probabilities 


Eq.  (9),  shown  below,  was  found  by  solving 
Eq.  (4)  for  Ki  and  substituting  relative 
frequencies  for  probabilities. 


Ki-  |log  [ RF (WaWb) +RF (W^) +RF (W^) -l]  -log 


{ 


1 

log  RF(iJjj) 


} 


(9) 


RF(Wa)  > RF(w|,) 

The  relative  frequencies  in  Table  3 and  the 
complements  of  the  values  found  in  Table  1, 


RF (S)-l-RF (W) , were  substituted  into  Eq.  (9)  and 
thirty-six  solutions  were  found  for  Kj,  one  for 
each  pair  of  sites. 


The  K}  values  are  given  in  Table  4 and 
plotted  in  Fig.  2.  A curve  of  the  form 
a-  D** 

K-  - — _ (10) 

P bec^  D* 

was  fitted  to  the  data  points. 

TafcU  4.  Solutions  (or  Kj  obtained  from  Eq.  (9)  in  winter. 
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uc 


• 164 
.226 


.398 

.378 

.371 


.558  .825 

.547  .417 

. '.39  .812 

.3t7  .455 

.250 


.822  .702 

.818  .895 

.809  .887 

.457  .574 

.235  .458 


.845 

.843 

.838 

.770 

.885 


.171 


.404 


Figure  2.  The  parameter  Ki  plotted  as  a function 
of  distance.  The  dots  and  X's  are  solutions  to 
Eq.  (9)  for  winter  and  summer,  respectively. 

The  curves  are  of  Eq.  (14)  with  a^-l . 2 in  winter 
(solid)  and  aj-2.7  in  suaner  (dashed). 

Eq.  (10)  was  chosen  because  of  its  charac- 
teristic shape  and  because  it  met  the  two 
necessary  conditions,  that:  (1)  equals  sero 

when  the  distance,  D,  equals  serot  and  (2)  by 
setting 

b - D*  (11) 

where  D is  the  distance  at  which  the  events  are 
statistically  independent,  Up  equals  one  when  D 
equals  08.  At  distance  Ds  when  equals  one 
Eq.  (10)  can  be  written  * 

> ■ > - (u) 

and  from  Eqs . (11)  and  (12)  it  can  be  seen  that 
cp  » % -1.  Thus,  if  the  distance  D,  is  known, 
a p ana  a are  the  only  parameters  that  require 
evaluation. 

The  parameter*  was  set  equal  to  one-half  and 
the  distance  Ds  was  estimated  to  be  650  and  550 
nautical  miles  in  winter  and  simmer,  respectively, 
based  on  extrapolation  of  the  points  in  Fig.  2. 

Eq.  (10)  was  solved  for  a^  using  the  above 
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information.  The  resulting  equation  is 
K-  (b-  f5) 


where  for  winter  b 
for  summer  b - J D~ 


The  36  Kj  values  given  in  Table  4 together 
with  the  distances  between  the  sites  were  used 
to  find  estimates  if  a^  for  winter.  The  average 
value  of  the  36  a^'s,  used  as  the  first  guess, 
was  1.19.  Because  the  points  are  unevenly  dis- 
tributed along  the  curve  a "best"  fit,  to  the 
authors'  eye,  was  found  by  slightly  varying  the 
a1  value  from  its  average.  The  value  chosen  was 
1.20.  The  susmer  average  was  2.69  and  the  value 
chosen  was  2.70. 


Substituting  Eqs.  (11)  and  (12)  into  (10) 


yields: 


3.5  Three-site  Estimated  Joint  Probabilities 


Eq.  (15) , shown  below,  was  found  by  solving 
Eq.  (6)  for  K]  and  substituting  relative 
frequencies  for  probabilities. 

K2-  {log  [ 1-RF (H^c) _RF ("*> -*F (^b) 

-RF  (5q)  +RF  (H^)  RF  (flg)  Kl+RF  (STa)  RF  ((£) Kl 

+RF  (WB)RF(j^,)KlJ  - log  RF  (5*)  RF  (1%)  KlJ- 


Table  5 gives  values  of  Ds  and  a 
winter  and  suseier.  1 


Curves  of  are  shown  in  Fig.  2 


l log  RF  (wc) J 
RF(5^)£  RF (5g)  ^ RF (W^) 


The  relative  frequencies  in  Table  6 (and  the 
others,  not  shown,  involving  sites  other  than 
LGA) , the  complenent8  of  the  values  found  in 
Table  1,  and  values  using  a^  equal  to  1.2  and 
2.7  were  substituted  in  Eq.  (15)  and  eighty-four 
solutions  were  found  for  *2,  in  winter  and  sumner 
one  for  each  set  of  three  sites. 


The  solutions  are  plotted  in  Fig.  3 and 
sample  solutions  are  shown  in  Table  7.  A good  : 
to  the  data  points  was  obtained  when  a2  was  set 
equal  to  1.00  and  2.30  in  winter  and  susmer, 
respectively. 


With  a ^ equal  to  1.20,  Eq.  (14)  was  substi' 
tuted  into  Eq.  (4)  and  estimated  joint 
probabilities  were  found  for  winter.  They  are 
shown  in  Table  3.  The  agreement 'between  the 
relative  frequencies  and  the  probabilities 
estimated  by  the  model  is  excellent. 


3.4  Three-site  Joint  Relative  Frequencies 


DISTANCE  (mtIM) 

Figure  3.  The  parameter  K2  plotted  as  a function 
of  distance.  The  dots  and  X's  are  solutions  to 
Eq.  (15)  for  winter  and  sunner,  respectively. 

The  curves  are  of  Eq.  (14)  with  a2-1.0  in  winter 
(solid)  and  a2=2. 3 in  stumer  (dashed). 


The  nine  sites  provided  84  combinations  of 
three-site  concurrent  precipitation  occurrences 
Sample  relative  frequencies  of  these  three-site 
joint  occurrences  are  shown  in  Table  6. 


- 

WL 

ML 

tot 

DCA 

UC 

■DU 

JTt 

.113 

,0M) 

.0871 

.0585 

.0585 

.0500 

.0322 

(.117) 

(.0869) 

< .0870) 

(.0620) 

(0616) 

(.0501) 

(.0322) 

wm 

.0892 

.0891 

.0802 

.0802 

.0517 

.0331 

(.0915 

(.0700) 

( .0845) 

(.0842) 

(.0515) 

(.0331) 

no. 

.0857 

.0573 

0589 

.0482 

.0298 

(.0825) 

(.0575) 

( .0589) 

( .0447) 

(.0288) 

ML 

.0581 

0590 

.0453 

.0271 

(.0588) 

(.0583) 

(.044:> 

(.0257) 

AON 

.0558 

.0438 

.0258 

(.0582) 

(.0444) 

(.0255) 

OCA 

.0438 

.0258 

(.0443) 

(.0258) 

uc 

.0273 

(.0289) 

Winter 

u 

! 

0 

650 

550 

*1 

1.20 

2.70 

•2 

X a 00 

2.30 

•3 

0.90 

2.15 

•4 

0.85 

2.10 

*5 

0.822 

2.08 

•6 

0.803 

2.068 

•7 

0.790 

2.062 

•s 

0.780 

2.057 

490 


T*bl*  7.  Mapl*  solution*  for  Ea  ofet*ln*d  turn  Eq.  (IS)  b***4  on 

rolotlvo  fr*qu*nci«*  of  joint  occurrence*  of  precipitation 
•t  USi  *nd  **ct»  of  tte  otter  toonty-oiqht  p*lr*  of 
ataorvatloo  sit**  in  winter. 


« 

ML 

BAL 

KM 

DCA 

RXC 

MNJ 

JM 

.122 

• 111 

• Hi 

.099 

.091 

.117 

.114 

— 

.139 

• 149 

.122 

.125 

.159 

.157 

ML 

.409 

.390 

.392 

.412 

.410 

IN. 

.530 

.542 

.542 

.545 

ACM 

.544 

.542 

.549 

MB 

.549 

.576 

I1C 

.47J 

and  K2  values  were  substituted  into 
Eq.  (6)  and  estimated  3-site  joint  probabilities 
were  found  for  winter.  A sample  of  them  is 
shown  in  Table  6.  There  is  good  agreement 
between  the  observed  and  estimated  values. 


3.6  Four-site  to  Nine-site  Joint  Relative 

Frequencies 

The  nine  sites  provided  126,  126,  84,  36,  9 
and  1 combinations  of  4,  5,  6,  7,  8 and  9-site 
concurrent  precipitation  occurrences,  respective- 
ly. Sample  relative  frequencies  of  these  four- 
site  to  nine-site  joint  occurrences  are  shown  in 
Table  8.  In  order  to  conserve  space  only  a few 
of  the  possible  382  joint  occurrences  are  shown. 
They  are  believed  to  be  typical  and  include  the 
highest  and  lowest  observed  relative  frequencies. 


Kbit  I.  1**1«  relative  frequencies  of  joint  occurrences  of  precipitation  at  foes  to 
si*  sites  and  corresponding  estimated  probabilities  obtained  froa  the  — Jr 
in  winter.  Xn  order  to  conserve  specs  only  a few  seta  of  aitee  are  included. 
They  include  a range  of  observed  relative  frequencies  fra*  the  lowest  to  the 

highest. 


Sites 


in  kv  m lu 

KXJ  LGA  m JFK 

AIM  RXC  LGA  JFK 

MM  KL  RXC  LG* 

MM  PNL  LGi  (W 

MM  LGA  B*  JFK 


■inter 

Observed  Calculated 


.02)4  .0214 
.0305  .0)04 
.0)99  .0422 
.0401  .0401 
.0546  .0590 
.0561  .0620 


mm  nu  me 
M mi  ML 
mm  mil  me 
AIM  no.  OCA 
MM  PNL  LGA 
PNL  LGA  EWR 


LGA  JFK 
me  LGA 
DCA  LGA 

LGA  EWR 


BAL  JTT 


.0200  .0221 

.0206  .0204 

.0)70  .0Mb 

.04*1  .0557 

.0510  .0561 

.05*1  .06)1 


■ou  PNL  me  OCA 

AIM  ROU  PNL  me 

AOW  RDO  PNL  LGA 

mm  pnl  me  DCA 

ACM  PNL  OCR  LGA 

PNL  DCA  LGA  CNR 


LGA  JFK 
lga  am 
on  JTK 
LGA  JTT 
RAL  JTT 
BAL  JFK 


.01*0  .0195 
.0190  .0214 
.0212  .0217 
.0)44  .0M1 
.0452  .0512 
.04*4  .0544 


3.7  Four-site  to  Nine-site  Estimated  Joint 
Probabilities 


Relative  frequencies  were  substituted  for 
the  probabilities  in  Eq.  (8)  and  solutions  were 
found  for  all  126  combinations  of  4-site  joint 
precipitation  occurrences.  The  solutions  are 
plotted  in  Fig.  4.  a good  fit  to  the  data  points 
was  obtained  when  a3  was  set  equal  to  0.90  and 
2.15  in  winter  and  summer,  respectively. 

Relative  frequencies  were  substituted  for 
the  probabilities  in  Eq.  (8)  and  solutions  were 
found  for  all  126  combinations  of  5-site  joint 
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Figure  4.  The  parameter  K3  plotted  as  a function 
of  distance.  The  lower  points  are  based  on 
winter  data  and  the  upper  points  on  summer  data. 
The  curves  are  of  Eq.  (14) . 

precipitation  occurrences.  The  solutions  are 
plotted  in  Fig.  5.  A good  fit  to  the  data  points 
was  obtained  when  84  was  set  equal  to  0.85  and 
2.10  in  winter  and  sianmer,  respectively. 


Figure  5.  The  parameter  K4  plotted  as  a function 
of  distance.  The  lower  points  are  based  on 
winter  data  and  the  upper  points  on  summer  data. 
The  curves  are  of  Eq.  (14) . 

The  curves  for  Itj,  K2,  K3  and  K4  for  both 
winter  and  simmer  are  shown  in  Fig.  6.  The 
separation  between  the  curves  decreases  with 
increasing  subscript.  In  other  words,  the  value 
of  "a",  the  only  parameter  permitted  to  vary  with 
the  number  of  observation  sites,  n,  decreases  less 
with  the  addition  of  each  new  site.  The  a^ 
through  84  values  are  plotted  in  Fig.  7. 


Figure  6.  Curves  of  I,,  Kj,  Kj  and  K4  for 
winter  (lower  set)  and  summer  (upper  set) . 

Curves  were  subjectively  drawn  through  the 
points  in  Fig.  7 and  extrapolated  to  a total  of 
nine  sites.  Values  for  aj  through  ag  were 
estimated  from  these  curves  to  obtain  values  for 
use  in  solving  for  K-,  Ke,  K7  end  Kg.  The  "a  - " 
values  are  shown  in  Table  5.  “ 
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Figure  7.  The  parameter  a plotted  as  a function 
of  the  number  of  observation  sites  included  in 
the  joint  probability  estimate.  The  upper  curve 
is  for  winter  and  the  lower  curve  for  sumner. 


Eq.  (8)  was  solved  to  obtain  joint  probabil' 
ity  estimates  for  all  4-  to  9-site  combinations 
of  sites.  The  examples  shown  in  Table  8 were 
selected  to  illustrate  the  errors  in  estimating 
the  lowest  and  highest  observed  relative 
frequencies  and  some  in  the  middle  of  the  range. 
There  was  generally  good  agreement  between  the 
observed  relative  frequencies  and  the  calculated 
estimated  probabilities. 
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Figure  8.  Curves  of  obtained  from  Fig.  2.  The 
dots  and  X's  are  solutions  to  Eq.  (9)  for  winter 
and  summer,  respectively,  when  recurrence 
relative  frequencies  are  substituted  for  joint 
occurrence  relative  frequencies.  A distance  of 
27  miles,  in  winter,  and  22  miles,  in  susmer,  was 
assumed  to  be  equivalent  to  one  hour  of  time  lag. 


At  the  conclusion  of  the  application  of  the 
model  to  spatial  joint  occurrences  a decision 
was  made  to  test  the  model  on  hourly  recurrence 
of  precipitation.  Only  sets  of  two  events,  at 
time  t and  t -t  x hours,  were  investigated  with 
the  following  results. 


30.  Again  there  is  good  agreement  between  the 
observed  relative  frequencies  and  the  calculated 
recurrence  probabilities. 


4.1  Recurrence  Relative  Frequencies 


Recurrence  relative  frequencies  of  precipi- 
tation for  the  nine  East -Coast  stations  are 
shown  in  Figures  published  by  Lund  and  Grantham 
(1976) . They  are  tabulated  for  selected  hours 
from  one  through  30  in  Table  9.  Eq.  (9)  was 
solved  for  using  the  relative  frequencies 
given  in  this  table.  It  was  found  that  correla- 
tion between  precipitation  occurrences  decays 
about  the  same  amount  in  one  hour  as  it  does  in 
a distance  of  22  miles,  in  summer,  and  27  miles, 
in  winter.  The  "K^"  curves  shown  in  Fig.  2 are 
repeated  in  Fig.  8 and  Kj  values  obtained  from 
the  recurrence  relative  frequencies  are  plotted 
in  Fig.  8.  The  scale  is  adjusted  so  that  22 
miles  equals  one  hour  in  summer,  and,  27  miles 
equals  one  hour  in  winter.  The  agreement  between 
spatial  Kj  values  and  temporal  K1  values  is 
excellent. 


The  model  developed  for  estimating  joint  prob- 
abilities of  weather  events  is  intended  for 
general  usage.  However,  the  parameters  of  the 
model  must  be  obtained  either  from  data  or  theory. 
They  will  vary  with  the  unconditional  probability 
of  the  event  and  the  spatial  decay  of  correlation. 
The  test  of  the  model  on  winter  and  summer  joint- 
occurrences  of  precipitation  indicates  that  the 
parameters  are  trail  behaved.  Because  this  is  a 
dependent  sample  test,  the  general  applicability 
of  the  model  parameters,  even  for  precipitation 
occurrences,  is  not  known. 


The  model  is  being  tested  on  other  weather 
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Eq.  (4)  was  solved  for  probabilities  every 
22  miles  in  summer  and  27  miles  in  winter.  The 
probability  estimates  are  given  in  Table  9 
corresponding  to  selected  hours  from  one  through 
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